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ABSTRACT 


Important insights into the dynamics of spherically symmetric AdS-scalar held per¬ 
turbations can be obtained by considering a simplihed time-averaged theory accurately 
describing perturbations of amplitude e on time-scales of order 1/e^. The coefficients of the 
time-averaged equations are complicated expressions in terms of the AdS scalar held mode 
functions, which are in turn related to the Jacobi polynomials. We analyze the behavior 
of these coefficients for high frequency modes. The resulting asymptotics can be useful for 
understanding the properties of the hnite-time singularity in solutions of the time-averaged 
theory recently reported in the literature. We highlight, in particular, the gauge depen¬ 
dence of this asymptotics, with respect to the two most commonly used gauges. The harsher 
growth of the coefficients at large frequencies in higher-dimensional AdS suggests strength¬ 
ening of turbulent instabilities in higher dimensions. In the course of our derivations, we 
arrive at recursive relations for the coefficients of the time-averaged theory that are likely 
to be useful for evaluating them more efficiently in numerical simulations. 



1 Introduction 


Anti-de Sitter (AdS) space-time is linearly stable. Yet, non-linear dynamics of its pertur¬ 
bations is characterized by complex interplay of stable and unstable behavior. Numerical 
simulations for a massless, spherically symmetric scalar held [1-6] indicate that, while some 
initial data of small amplitude e collapse, after multiple rehections from the AdS boundary, 
to form black holes on time scales of order 1/e^, other small amplitude initial data remain 
small for much longer (possibly, for all) times; see also [7-9] for generalizations to massive 
scalar helds. 

Attempts to control the small amplitude dynamics in terms of naive perturbative expan¬ 
sions in the perturbation amplitude £ are plagued by a profusion of secular terms [1] growing 
in time and invalidating the perturbative expansions. This profusion is due to the perfectly 
commensurate spectrum of frequencies displayed by the linearized AdS perturbations. The 
perturbation theory can be re-structured using a variety of techniques equivalent at the 
lowest order (multiscale or renormalization group resummation, averaging), producing an 
improved perturbative expansion accurately describing the small amplitude dynamics on 
time scales of order [10-14]. 

As a result of applying the resummation techniques in the small amplitude regime, one 
obtains a system of effective equations governing slow variations of the amplitudes and 
phases of the AdS normal modes under the effect of non-linearities. We usually refer to 
these equations as “(renormalization) how equations” or “time-averaged equations”; they 
have also been called “two-time framework equations” [10] and “resonant system” [13]. The 
coefficients in this system are complicated expressions in terms of integrals of the AdS mode 
functions and their derivatives, obtained in [11,12]. 

In [13], the how equations were solved numerically for AdSs. It was observed that, for a 
variety of initial data for which collapse occurs in the fully non-linear description, a singu¬ 
larity in solutions develops in hnite time. At this singular moment, the time derivatives of 
the phases diverge logarithmically, while the amplitude spectrum, initially rapidly decreas¬ 
ing in the ultraviolet, attains a power-law form 1/n^, where n is the mode number. This 
signihes efficient weakly turbulent transfer of energy to short wavelength modes, which is 
the main mechanism behind the collapse dynamics. 

Since the key questions on the fate of AdS perturbations involve ultraviolet dynamics, 
it is important to develop asymptotic techniques for analyzing the high-frequency contribu¬ 
tions to the flow equations. Here, we initiate such analysis by addressing the high-frequency 
asymptotics of the interaction coefficients appearing in the flow equations (see also our par¬ 
allel work [15]).^ To this end, we hrst review, in section 2, the basics of constructing the 
flow equations for spherically symmetric AdS-scalar held perturbations. In section 3, we 
develop iterative relations for the interaction coefficients in these equations. These iterative 
relations may be of independent value for the purposes of optimizing numerical simulations 
of the how equations, as in [13]. In section 4, we use the iterative relations to reveal the 
power law behavior of the interation coefficients for large mode numbers, and emphasize 

^Some further considerations of ultraviolet properties of AdS perturbations were given in [18] after the 
present article had appeared as a preprint. 
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that this asymptotics is substantially affected by the time parametrization conventions (two 
gauges, corresponding to labelling the time slices by the proper time in the center and at 
the boundary of the perturbed AdS, have commonly appeared in the literature). We con¬ 
clude with a discussion of the qualitative implications of the power law behavior we find, in 
particular, the qualitative dependence of the near-AdS dynamics on the number of spatial 
dimensions. 


2 Weakly non-linear AdS dynamics and the effective time-averaged 
description 


Reviews of the general spherically symmetric AdS-scalar field set-up are widely available in 
the literature, for example [2] or our own work [11,12]. We shall therefore limit ourselves 
to a very basic summary of the relevant formulae. 

One considers the following spherically symmetric metric: 

^ ^ - Ae~'^^dt^ + sin^ x , (1) 

cos^x \ A I 


with A{x,t), 6{x,t) and a scalar field (f>{x,t) depending on the time coordinate t and the 
radial coordinate x G [0,7r/2). A = l, d = 0 = O corresponds to empty AdS^+i. 6{x,t) 
has some residual gauge freedom 6{x,t) t S{x,t) + Qit), corresponding to redefinition of 
time. Two possible gauge fixing conditions have appeared in the literature: d(0,t) = 0 
[1, 10, 11,13, 14] and 5{7i/2,t) = 0 [5, 12]. We shall refer to them as the ‘interior’ and 
‘boundary’ time, respectively. 

The equations of motion take the form 


$ = (Ae-‘^n)', 

A' = ^ (A - 1) - /rz/ (<l>2 + n^) A, 



(2a) 

(4= + n=), 

(2b) 


A = -2fiuA^e-^<l>U, (2c) 

with <h = 0' and fl = A~^e^(j) (where dots and primes denote the t- and ^-derivatives, 
respectively). 


jji[x) = (tanx)”* ^ 


and 


p^x) 


{d-l) 

/i'(x) 


sin X cos X 
(tana;)‘^“^ 


(3) 


These equations, as they are, are forbiddingly complicated. Further analytic progress can 
be made by considering small amplitude solutions. Naive power law expansions in the field 
amplitude result, however, in unacceptable perturbation series plagued by growing ‘secular’ 
terms [1,10-12]. Improved perturbative expansions can be constructed using multiscale [10], 
renormalization [11] and averaging [12] methods, all of which are equivalent at the lowest 
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non-trivial order. The general idea is to use as the starting point for weak held expansions 
a linearized solution with slowly modulated integration constants, namely 


OO 

0(x, t) = An{eH) en(x) cos(a;„t + Bn{eH)). (4) 

n=0 

Here, e„ are the mode functions of a linearized massless scalar held in hxed AdS^+i, 


(X 


kn 


COS X 


9 ^■’9. J 


(cos(2x)) with 


kn — 


2 y/r?T(nr+”d^njI 

r (« + i) 


( 8 ) 


P!f'’^\x) is a Jacobi polynomial of order n. Un = d + 2n is the famously integer, fully com¬ 
mensurate spectum of the linearized normal mode frequencies in AdS, largely responsible 
for the rich non-linear dynamics of AdS perturbations. 

The ehective equations for the slow time dependences of An{eH) and Bn{eH) can be 
deduced using a number of equivalent methods [10-12]. Their analytic form was presented 
in [11,12]. The equations are known to describe the exact non-linear dynamics accurately 
on time scales of order 1/e^, if the scalar held amplitude is of order £ [12]. Some phe¬ 
nomenological consequences of these equations have been explored in [13,14]. 

The averaged equations [12] are of the form 


2uji dAi 
dt 


2uiAi dBi 

dt 


{hj} {k,l} 

- ^ ^ ^ SijkiAiAjAk sm{Bi + B^- Bi- Bj), 

^ i j k ^ 

idi -\-(di 

i^l {i,j} ^ {k,l} 

— TiAf — ^ RiiA^Ai — ^ ^ ^ SijkiAiAjAk cos{Bi + B^ — Bi — Bj 

i i j k 

~\~^j 


( 6 ) 


(7) 


Note that a priori resonant frequency quartets other than Ui + Uj = Uk + ooi might have 
appeared in the sum, but they do not in fact contribute [11] due to special selection rules 
in AdS (see also [16,17]). The coefficients of the terms that appear in the time-averaged 
equations can be expressed in terms of certain integrals of products of the mode functions. 

In boundary time gauge (J(7r/2) = 0), in which the slow dynamics turns out to be 
manifestly Hamiltonian [12], the coefficients take the following form: 


T; = 


nil + ^Yiiii -F 2ui 


( 0 , 0 ) 


2uiWi 


( 1 , 0 ) 

nil y 


( 8 ) 


Rii =- 


uf + uf 






{ufX., 


illi 




un) + 2 




11 nil 


ojtYim 


UJt 


OJJ 


+ '^{Yiui J- Yiiii) 


UJt 


■ UJ 


(A',„ - A'„„) + 


( 0 , 0 )^ 


41,0) 


41,0) 
all ’ 


(9) 
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Sijkl 


1 1 
+ 


+ 


4 \Ui + Uj Ui — uik ujj — ujk 


1 1 
+ 


4 \u}i + ujj Ui — Uk Uj — Uk 

1 1 

+ 


4 \Ui + Uj Uj — Uk Uj — Uk 
1 


4 \Ui + Uj Ui — Uk Uj — Uk 
The integrals that appear in these expressions are dehned by 

^7r/2 


iy^i^k^lXjikl ^j^ijkl) 

j(^lXkijl ^l^ikjl)' 


( 10 ) 


Xijkl 

^ijkl 






r‘7r/2 


t/2 


r‘7r/2 


dx e'i{x)ej{x)ek{x)ei{x){^i{x)fv{x), 
dx ei{x)ej{x)e'k{x)e'i{x){^i{x))^iy{x), 

I pX 

dxei(x)ej(x)/i(x)i/(x) / dyek{y)ei{y)ii{y), 

Jo 

I nx 

<\xe[(x)e'Ax)iJ,(x)v(x) / Ayei(y)ei(y)lx(y). 


(lla) 

(llb) 

(llc) 
(lid) 


Note that the coefficient Rij = Rji is symmetric. Whenever the resonance condition Ui + 
Uj = Uk+ ui is satished, one also has Sijki = Sjiki and Sijki = Skiij- 

In the interior time gange (5(0, t) = 0), the expressions for R and Rij are slightly 
more complicated [11,12], with the additional terms uf{Aii + u^Vu) and u'^i^Aa + u'^Vu), 
respectively. The V and A coefficients are dehned as. 


f7r/2 




dx ej(x)ej(x)p(x)z/(x). 


-71/2 


Ajj 


dx e'(x)e((x)p(x)z/(x). 


( 12 ) 


Note that the difference between the two ganges only affects the eqnation for the phases 
Bn, bnt not the eqnations for the amplitndes An- In this gange, the R coefficients are not 
symmetric and the slow dynamics is not manifestly Hamiltonian [12]. 


3 Recursive analysis of the interaction coefficients 

3.1 Recurrence relations for the mode functions 

We hrst analyze the mode-nnmber strnctnre of the AdS^+i mode fnnctions given by (5). 
Using the recnrrence relations for the Jacobi polynomials Pi“’^^(x), 

2(?t, + l){n + a + /3 + l)(2n + a + /9)P,^“f^(x) = —2{n + a){n + /3){2n + a + f3 + 2)P,^"f^(x) 
+ (2n + a + /I + 1) |(2n + a + /I + 2)(2?7. + a + /3)x + o? — P^’^’^^(x) 

(13) 
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and the derivative formula, 


(2n + a + /3 + 2)(1 - = -2{n + l){n + a + (3 + 

P {n -\- Oi -\- j3 -\-V) (ex — /5 + (2?r + a + /3 + 2)x) (^x) (14) 

we can show that the mode functions satisfy the following identities, 

fi{x)P{x)en{x) = -{d - 1 ) ( 1 -y—r ] (^n{x) + A/(n + l)(n + d) ( —^ ) en+i{x) 

V ^n-1/ \Un + lJ 

+ ^/n{n + d-l) ( — ^—-]en-i{x) (15) 

\UJn ^ / 


and 

2ii{x)v{x)Er,{x) = -{d- 1) x ) A/(n + l)(n + d) ^ en+i{x) 

- ^/n{n + d-l) (—^ ) e„_i(a;). (16) 

\0Jn -L / 

We have introduced the notation En{x) = {l/ujn)e'^(yx). One can show that En{x) is also 
an orthonormal family of functions with weight fi{x). By taking derivatives of these and 
using —{jxe'^y = uj'^fien and — (/iz/')' = 4/iz/, we can obtain the relations 

n{x)P{x)En{x) = - (d - 1) ( 1 + d 1 1 En{x) + \/{u + l){n + d) ( — \-r] En+i{x) 

V ^n-1/ \Un + lJ 

+ ^/n{n + d-l) ( — En-i{x) (17) 

\UJn -L / 


and 


2n(x)u(x)en(x) = - (d - 1) E^i^x) - ^/{n + l){n + d) ^ En+i{x) 

+ ^/n{n + d - 1) (—^ ) P„_i(a;). (18) 

\Un - 1/ 

Note that up to a few signs, these recurrence relations preserve their form under interchange 
of en{x) and En{x). This feature will allow us to streamline some of our derivations below. 

3.2 Recurrence relations for the X coefficients 

A useful simplification is attained by expressing the coefficients 

TT 

Xijki = / dix{^i{x)v{x)e({x))ej{x)ek{x)ei{x)^{x) (19) 

Jo 


5 



in terms of the totally symmetric combinations 


Xijkl 


dx ei(x)ej(x)ek(x)ei(x)jj,(x) 


( 20 ) 


with the help of the recnrrence relation (16). This resnlts in the formnla 


T \ f 

—-^nmpq — J-j I 

OJn ^ \ 


\nmpq “t~ 2 '\/(^ “1“ l)(^ ^ 

- ]^^n{n + d-l) 


a;„ + 1 

1 


X{n+l)mpq 
X(n—l)mpq- (21) 


We shall now derive a set of recnrsive relations for x. Once y have been determined, X 
can be extracted nsing the above formula. 

Consider hrst the product fiu'enem- One can apply the recurrence relation (15) either 
to ^i^'cn or to resulting in two different re-writings of this product. Equating these 

two different re-writings, multiplying by ^CpCq and integrating over x, one obtains 


{d - 1) Xnmpq + {m + l){m + d) + 7^771(771 -Fd- 1) 

= (d-l) 


X.n{m—l)pq 
^m, 1 


,W„+\/(" + l)(n + <i)(^^^^) + yn(n + <i-l) (*5^) , 

( 22 ) 


(By total symmetry of y, a similar relation holds for any other pair of indices.) 
Furthermore, one can consider the relation 



dx jX V C^C^CpCq 


-2(d-l)y, 


nmpq 


X. 


nmpq 


X, 


mnpq 


X, 


pnmq 


X, 


qnmp' 


(23) 


Here we have used integration by parts moving the derivative off z/', and the fact that 
/iV = (d — 1). This relation can be easily converted in a closed iterative relation for y 
alone. On the right-hand side, we can eliminate X in favor of y with (21). On the left-hand 
side, applying (15) to pz/'cn re-expresses everything through y alone. The result is 


+ 1 


X{n+l)mpq 


(d - 1 ) ~ \ V{n + ^){n + d) 

-f ^3/77(77 + d-l) ^ X(n-l)mpg| + {n ^ TTl} + {u-H- p} + {n ^ q} 

(d ^) (^ 2 ^ Xnmpq T \/('^ T 1) (77 -|- d) f j —^ X{n+l)mpq 


^/n{n + d-l) ( —) 

\U!n 1 / 


X{n—l)mpq- 


(24) 


We shall now present the analogous relations for the H-coefhcients, and then explain how 
(22) and (24), and the corresponding equations pertinent to the X-coefficients, can be used 
to set up effective recursive evaluation algorithms. 
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3.3 Recurrence relations for the Y coefficients 


We have already remarked on the immediate parallels between the recurrence relations 
for en{x) and En{x) = (l/a;„)e^(x). These parallels make the analysis of the Y coefficients 
essentially identical to what we have presented above for the X coefficients. Using recursion 
formula (18), we can write 


■y, _ 

^ ijkl 


da: el{x)ej{x)e'y)e'i(x)(fi(x)fv(x) 


(25) 


in terms of 

as follows: 

^nmpq 

(jJn^p^q 2 


Ajki= / dx Ei{x)Ej{x)Ek{x)Ei{x)iJ,{x) 


Ur. 


2 j '^nmpq „ (jTl T l)('m. T ctj 

OJxrr, -L / ^ 


+ 1 

1 


(26) 


V^n(m+l)pg 


+ + d-l) ( -- 1 'lpn{m-l)pq 

^ \Um. — i 


(27) 


Considering the action of (17) on ^u'EnEm results in 
1 


-(d-1) 

= -id-l) 




cut 


+ 1)(" + d) (+ y„(„ + d - 1) ( W')""' 
V “T -L / \ 

+ y(m+l)(m + d) ( + v/m(m + d-l) 

\ “r i / 


U)m 1 


(28) 


On the other hand, one obtains by integration by parts 


dx fi^u'E^Err^EpE, = 2{d- l)i;r^mpq + + ^mYnmpq ^ UJpY„ ^ UqY^ 


'0 


^n^p^q ^n^mP^q ^n^rnP^p 


Using (17) and (27), we can conclude that 

1 


-Pd-l) 




+ ^y/n(n + d - 1) ^ 

= -(d-l) 


'4^nmpq „ A^(n + 1) (?7. + (7) i i 1 ) ^( n+l)mpq 
^ \^n + J- / 

+ {n ^ m} + {n ^ k} + {n ^ 1} 




+ 3^ V’ nmpq + y(»+l)(» + :i) (;^) V', n+l)mpq 


+ y/n(nTd^^ ( —- ) ?/>(„ 

— l)mpq 

\UJn “ / 


(30) 
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3.4 Recursive evaluation of the X and Y coefficients 


One important feature of the recursive expressions in the previous section, is that we can 
easily compute all the X and Y coeffients recursively. We will explain this in detail for the 
X coefficients. The procedure is completely analogous for the Y coefficients. 

One starts by evaluating y as defined by (20). For Xnmpq, we call L = n + m+ p + q 
the ‘level’ of that coefficient. At level 0, there is only one coefficient, which can easily be 
evaluated analytically. 


Xoooo — 


(r(d))^ r(f) 
r{2d) f 


(31) 


Now note that the formula (24) relates coefficients of level L—1 and L to multiple coefficients 
of level L + 1. However, by substituting equation (22) three times in (24), we can find a 
formula that expresses a single coefficient of level L + 1 in terms of coefficients of level L 
and L — 1: 


An I iV-n I ^^ i^n + UJm+(^p + (^q + ‘^) 

-y y 7 l l)yn d) , , -I \ X(n-\-l)mpq 

[UJn -r 1 ) 


= {d-l) 


tat 


2 2 

H-H-— iXm + cUp + cUg + 1) - 




OJm. 1 1 


- 1 ) 


Xnmpq 


+ y/n(n + d- 1)—— ^ —^X(n-i)mp 5 + \/m(m + d - 1) ( ^ 

[pJri t J y (Urn t J 

+ y/p(p + d- 1) Xnm(p-I)q + ^/qiq + d - 1) 


Xnmp(q—1) ■ 


Xn(m—l)pq 


(32) 


So once all coefficients Xnmpq of two adjacent levels are known, one can easily compute all 
coefficients of the next level. Since the need for explicit integration is completely eliminated, 
no fast oscillating integrals have to be computed explicitly, and this method will be more 
efficient than a direct evaluation of the integrals of mode functions. After the recursive 
computation of the y coefficients, one only needs to use (21) to determine the X coefficients. 

The method for the Y coefficients and the p coefficients is completely analogous. One 
simply needs to combine equations (28) and(30) and start with the value^ 


Aoooo 


8(r(rf))^ r(f 

r(2(i+1) 


-i)r(i + 2) 

(r(i)r 


(33) 


Once all X and Y coefficients are known, we can determine the coefficients Sijki that 
appear in the flow equation (6) for the amplitudes. In order to determine the Rij and 
Tj coefficients, we also need to compute the Wiju integrals appearing in (11c) and (Hd). 
Though there are much fewer R and T coefficients than S coefficients, and hence optimizing 
their evaluation is probably less crucial, we shall now explain briefly how they can be 
efficiently computed as well. 

^There were typos in this formula in the original version of the paper. We thank Brad Cownden and 
Andrew Frey for bringing that to our attention. 



















3.5 Computation of the W coefficients 

Consider the PC-coefficients as defined in [12], 


= 


r-7r/2 




In appendix D of [12] we established the relations 


( 34 ) 


= W.,,, 

H - ^f) < 1 ?' = Xuj, - 

H - 



Ol - 

kijl 1 

{u,l - a,?) <«' 

^kijl 7 

{iol - 0,?) <):;' 


\Y 

^ iljki 

(35a) 

^ikjl ^iljk^ 

(35b) 


(35c) 


These allow us to express all for which fc 7 ^ / in terms of X and Y coefficients. 

However, in the expressions for Rij and R we need the K-coefficients with k = 1. Assuming 
that we already know the y, 'ip, X and Y coefficients, these can be calculated recursively 
in the following way. 

Consider the integral dx ejCj/iz/ dy fi{fiu')ekei. By letting the (/in') work on e*, 
and ei respectively and using the recurrence relations (16), we find that 




+ ^/k{k + d-l)( 

\ Uk — 1 J 

(36) 


Now we substitute / = A; +1 in this formula and replace all wi^Jq with p ^ q that appear in 
the remaining expression by (35). What remains is a recursive relation that relates W-j^^ 
simply to This allows for a fast recursive computation of the coefficients, 

once the value of is known. 

The most efficient way to compute 1T/^° for all i and j is also by a recursive procedure. 
In close analogy to the other derivations we have presented, one obtains 


(d-1) 


a;/-l 


+ \/(* + 1 )(^ + d) 


W, 


( 0 , 0 ) 




(i+l)jkl \ 


+ 1 / 


+ \/i(i + d - 1) 


W, 


( 0 , 0 ) 


(i-l)jkl 


Ui - 1 


/ 1 \ _ \ _ 


-l)kl 

r 


(37) 
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and 




+ d 






I ~ 2 V^(^ + 1 )(* + 


Ui 


I {i+l)jkl 


■ 1 ) 


Ui 






^ \^l ^ / ) 

(<i-i)<'!;?’ + 2 (<i-i)(;^) 

2Vi(i + <i-l) (^) H 


+ 1 / 

<H- j} + {i -H- k} + {i 1} 

W'SS' + i^i.^+m + d) M/^y 


Using these relations, the W^jqq^ coefficients can be compnted in exactly the same way as 
for the X and i/j coefficients. Only now we take L = i + j to be the ‘level’ of 

We can easily extract the coefficients once all the are known. On the one 

hand, we consider with e, expressed as —{fxe^y/u^, and nse integration by parts to 

throw one of the two thns inserted derivatives off e*, obtaining the relation 

- - Xii„ (39) 

where we have dehned the integral dx e^ejfiu' dy fi. On the other 

hand, we have by integration by parts, 

+ USi = W.‘“' + 2(d - + A'yt, + AyH + Xui, + Xujf (40) 

Here we have nsed the fact that {yi^'Y = —4pz/ and the expression for the integral (23). 
Combining these two, we hnd that 

+ u] - = 2{d - l)x^Jkl + 2X,jki + 2Xj,ki + X^iji + Xu,k (41) 

Assnming all y, X and are known, this allows ns to compnte all 


3.6 Additional coefficients in interior time gauge 


If one works in the interior time (5(0, f) = 0), one also needs the Vij and Aij coefficients 
defined by (12). These coefficients can also be compnted by a simple recnrsive procednre. 

Jnst as before, we first consider the action of the recnrsion relation (15) either on yv'cn 
or on yv'cm in the prodnct yv'cnem- We then eqnate the results, multiply by yv and 
integrate to obtain 




(n—l)m 




UJr 


^^Vnm + \/{n + l)(n + d) + \/n{n + d - 1) ^ 

^^Vnm + V(m + l)(m + d) xjm{m + d-l)( 

V + J- / V 


1 

hn(m— 1 ) 


UJt 


UJri 


(42) 
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We then note that 


{d — l)Vij = {d — 1) / dixeiCj^v = / dx eiejfii'{^'u) = — dx fi 


dx {fiue'j)ej^Ty — / dx — / dx ejej/iz/(/iV) — 2 / dx (/iz/'ej)ej/ii/ 

Jo Jo Jo 

(43) 


Using the recnrrence relations resnlts in 

1 
2 


+ 1 


(d - 1) ^ V(" + l)(n + <i) ( 

+ 1 + rf - 1) ^ j V(„_i)m| + {« 44 m} 

= 2((i — 1) [ —^ j Vnm + 2 ^/{ji + l)(n + d) 


V( 


(nH-l)m 


I 


+ 1 


Ur 


(n+l)m 


+ 2'\/n(?7, + d — 1) ( -- J V(n-I)r 

\ t y 


(44) 


The above relations can be nsed to efficiently compnte the Vij coefficients, starting with 

V5,„ =-Hhi-(45) 

(<i + i)(r(i))^ 

In appendix B of [12], we established that 


= with Ci = 

Therefore, we can immediately obtain all ^-coefficients as well. 


_ 2y/d^ {i + d-l)\ 


ii 


(46) 


4 Ultraviolet asymptotics of the interaction coefficients 

We can extract the behavior of the interaction coefficients for large mode numbers by solving 
nrrence relations asymptotically in this limit. In particular, for the Xijki coefficients, 
.sider the recurrence relations (22) and (24). In the limit of large values of n, m, p 
these reduce to 


the recurrence relations asymptotically in 
we consider the recurrence relations (22) 
and q, these reduce to 

1 1 1 1 

\d IjXnmpg T 2^Xn{m+l)pq T ^Xn{m—l)pq \d IjXnmpg T 2 T(n+l)mpg T ^ 

and 

l)Xnmpq '^'^X(n+l)mpq T 2^T(n—l)mpij^ T ('^ 'kTlj {jl •<->■ p) -|- {jl ■<->■ 5) 

1 1 
\d IjXnmpg T 2 T(n-|-l)mpg T 2 


(48 
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By introducing finite difference operators Ak\nmpq = {X{n+i)mpq - X{n-i)mpq)/‘^ and 


^nXnmpq = X{n+i)mpq + X{n-i)mpq “ 2xnmpg , we cau Write these equatious as 

^mXnmpq ^nXnrapq 


(49) 


and 


nmpq + mA^2x nmpq +pA^ph nmpq +q^^qh nmpq l)Xnmpg 2 ^{,X.{n-\-l)mpq^X.{n—l)mpq\ 


(50) 

Note that our finite difference operators are simply proportional to discretized first and 
second derivatives. (49) is a discretized two-dimensional wave equation and can be solved 
by an ansatz of the following form (which takes into account the total symmetry of Xnmpq)'- 


Xnmpq = fiiu + ui + p + q) + / 2 (-n + ui + p + q) + f 2 {n -m + p + q) 

+ f 2 {n + m - p + q) + f 2 {n + m + p - q) 

+ fsin + m - p - q) + fs{n - m + p - q) + fsin - m - p + q), (51) 


where the function /a must be even, fs^—x) = fz{x). One can develop some intuition about 
(50) by likewise considering its continuum limit, assuming \xn+i — Xn\ -C \Xn\i which yields 
the Euler equation for homogeneous functions. 


ndnX + rndmX + pdpX + qdgX = {d - 2)x, (52) 

solved by taking fi{x) = AiX^~'^. (For d = 2, we expect a logarithmic behavior fi{x) = 
Ai\nx.) One can directly verify that such a solution is consistent with the original hnite 
difference equation (50) at large values of the indices. Overall, 

Xnmpq ^ A^ {{h + TTi - p - qY~‘^ + {ji - TU + p - qY~‘^ + {n - ui - p + qY~‘^) 

+ A 2 {{—n + m + p + qY~‘^ + {n — m + p + qY~‘^ + {n + m — p + qY~^ 

+ {n + m + p — qY~^) + Ai{n + m + p + qY~^- (53) 


(We have compared this expression to some explicit brute force calculations for special 
index values using Mathematica.) 

One can check that the asymptotic equations for the V'nmpg coefficients are exactly the 
same, thus we expect the same kind of behavior. Extending this to the X and Y coefficients, 
plugging this back in (10) and assuming that there are no non-trivial cancellations, we 
expect the following scaling 

Sxi,xj,xk,^i ~ ^ ( 54 ) 

for any i,j, k, I, with A 1. 

Similar considerations yield for the lE-coefficients 


W- 


( 0 , 0 ) 


M,Aj,Xk,Al 


A 


d-4 


W- 


( 1 , 0 ) 




Ad-2 


(55) 
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Hence, the R and T-coefficients in the boundary time gauge behave as 

Rxi,\j ~ Txi ~ A'^. (56) 

Note that the R and T coefficients in the interior time gauge feature extra terms, as 
explained above (12). Repeating the asymptotic analysis for the Vij coefficients, we hnd 
the equations 

n/^]^Vnm + md!^Vnm = {d - 3)Vnm (57) 

and 

(58) 

which are solved by Vnm ~ z^i(n + mY~^ + 102(^1 — mY~^. For d = 3, this goes to ln(n + m) 
and ln(n — m) instead, and for <7 = 2, this goes to a constant (the explicit value of this 
constant is given in appendix B of [12]). The dominant contribution to the T and R- 
coefficients becomes ojfVn and ufu’jVa, respectively, rather than coming from the terms 
already accounted for in (56), which become subleading in the interior gauge. For AdS 4 
this results in the behavior Tj In j and Rjn ~ In j, consistent with the one reported 
in [14]. For AdS^+i with d > 3, it results in the behavior Tj ~ and Rjn ~ 
consistent with the one reported in [13] from numerical evaluations in AdSs. Note that 
we do not hnd evidence for logarithmic modulations in the leading power-law behavior of 
the interior gauge interaction coefficients for d > 3 (a possibility suggested in footnote 9 
of [14]). 

5 Discussion 

We have presented some analytic considerations of the interaction coefficients appearing in 
the time-averaged equations describing slow energy transfer between the normal modes of 
a spherically symmetric AdS-scalar held system due to gravitational non-linearities. In the 
course of our considerations, we have developed an iterative procedure for evaluating the 
interaction coefficients, which is likely to present signihcant advantages over a brute force 
approach. Indeed, straightforward computation of the interaction coefficients can be done 
in essentially two ways, both of which are resource-intensive. Either one multiplies high- 
degree Jacobi polynomials symbolically and integrates the result exactly using computer 
algebra, or one resorts to a numerical treatment of the integrals of products of Jacobi 
polynomials, which become very rapidly oscillating as the mode number increases. Our 
recursive procedure eliminates the need for integration altogether (the only drawback is 
that one has to evaluate the interaction coefficients for all modes, rather than only for those 
that satisfy the resonance condition). This sort of optimization is likely to enable signihcant 
improvements in the numerical analysis of the time-averaged theory of the sort presented 
in [13]. It would also be interesting to know how the iterative evaluation competes in terms 
of computational efficiency with the explicit expressions for the interaction coefficients in 
the particular case of AdS 4 presented in [14], as those expressions are extremely long and 
require many hoating point operations per coefficient. 
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Our iterative formulas allow us to deduce the asymptotic power law dependences of the 
interaction coefficients. In the boundary time gauge, if all the mode numbers are scaled up 
uniformly controlled by a parameter A, all interaction coefficients scale as in AdS^+i. The 
ultraviolet growth of mode couplings becomes more harsh in higher dimensions, conforming 
to the strengthening of turbulent behavior in higher dimensions alluded to in [14], 

It would be interesting to explore how the power law (54) appearing in the scaling of the 

5 coefficients relates to the properties of the singular spectra developing in hnite time in 
solutions of the time-averaged system according to [13]. (This singular behavior is believed 
to reflect the black hole formation in the full non-linear theory, though establishing this 
relation rigorously would require going beyond the time-averaged description.) In [13], it 
was observed numerically that, in AdSs, the following behavior of the amplitude spectrum 
An emerges near a certain hnite ‘singular’ moment t*: 

An{eH) - (59) 

This solution was verihed to be crudely analytically consistent with the time-averaged 
amplitude equation (6) in the sense that both sides appear approximately independent of 
n at large n, given the scaling Sxi,\j,xk,\i ~ in AdSs. We hope that our results, when 
combined with generalizations of the numerical analysis of [13] to higher dimensions, will 
lead to a more detailed and more general understanding of such singular spectra. 

Another comment is that in [13], which used interior time gauge, the how equations for 
d = 4 led to hnite-time singular behavior of (the time-derivative of) the phases i?„, and this 
singular behavior seemed to be driven by the Rjn ~ rAj^ scaling. Our analysis shows that 
this scaling is due to a contribution that is absent in boundary time gauge, so it would be 
interesting to verify whether the singular behavior of the phases persists in boundary time 
gauge. 
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